Equivalence of the velocity and length gauge perturbation series 
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We derive a "master" perturbation expansion for the quantum transition amplitude in a light 
field between the field-free initial and final atomic states in the minimal-coupling (MC) "velocity" 
gauge. The result is used to prove that the traditional "velocity" and "length" gauge perturbation 
series are equivalent infinite series representations or branches of the same amplitude function, that 
are equal but in a common domain of convergence (if it exists). More generally, we show that they 
constitute only two members of a one-parameter family of infinitely many branches of the given 
transition amplitude. 
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A long-standing unresolved problem of the quantum 
mechanical perturbation theory of light-matter interac- 
tion is the existence of two distinct infinite perturba- 
tion series for a given transition amplitude in the so- 
called "velocity" and "length" gauges, that have hitherto 
resisted a demonstration of their mathematical equiva- 
lence. In addition, numerical calculations of the transi- 
tion probability based on the two series (albeit, for prac- 
tical reasons, only in their truncated forms) have fre- 
quently shown a significant discrepancy between them 
and/or with various experimental data. These and re- 
lated difficulties have led some authors to argue in favor 
of the length gauge (e.g. as opposed to the velocity 

gauge, since the former is based manifestly on a physi- 
cally "true" energy operator Q. Nevertheless, theprinci- 
ple of gauge invariance in quantum theory (e.g. 
requires that they aught to be equivalent. 

The purpose of this Letter is to derive a "master" per- 
turbation expansion of the quantum mechanical transi- 
tion amplitude in the minimal-coupling (MC) "velocity" 
gauge, and to use the result to demonstrate that the two 
perturbation series, that are traditionally obtained in the 
"velocity" and "length" gauges, are two branches of the 
same amplitude function and, hence, they are equal but 
in a common domain of convergence (if there is any). 
More generally, our result shows that they belong to a 
one-parameter family of infinitely many equivalent series 
representations (or branches) of the transition amplitude 
for a given transition process. 

The Schrodinger equation of an atomic system inter- 
acting with an electromagnetic field, in the minimal- 
coupling (MC) transverse gauge is given by: 



(in^-HMc{i))-^Mc{t) = 
ot 



(1) 



where, the total Hamiltonian of the interacting system is 



Hmc {t) 



2m 



+ eAQ{r) 



(2) 



with the four-potential = {Ao{r), A{t)) where, the 
scalar potential can be used to define the "atomic" po- 
tential, eAo{r) = Va{r), and A{t) is the transverse vec- 
tor potential of the light field. In the usual electric dipole 



approximation (e.g. Bohr-radius/ wavelength << 1), the 
vector potential A{t) depends only on t. As usual, the 
"atomic" Hamiltonian 



provides a complete set of eigen-states 



>< 



(3) 



(4) 



allj 



and eigen-energies, E^, which satisfy the eigenvalue equa- 
tion: 



>=i?j|0? >,aUj 



(5) 



We define the "atomic" Green's function Ga{t,t') by the 
equation 



Its solution is given by, 
n 



(6) 



^9it - t') > e-i^^(*-*') < </.j"| 

allj 



(7) 



This can be easily verified by its substitution in Eq. ^ 
and noting that the derivative of the theta-function is the 
delta function. Finally, we define, for later use, the total 
Green's function (or propagator) GMc{t,t') associated 
with the minimal-coupling Hamiltonian HMc{t)i by the 
inhomogeneous differential equation: 

{^iU^ _ HMc{t))GMcM = S{t - t'). (8) 
ot 

The total Hamiltonian HMc{i) can always be written 
as a sum of two terms, in infinitely many ways: 

HMcii) = Hs{t) + Vs{t);s = 1, 2, 3, • • • oo (9) 
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where, Vs(t) = HMc{t) — Hs{t), and the basis Hamil- 
tonians Hs{t) can be used to define the associated basis 
Green's functions, Gs{t,t'), by the equations: 

{in^^~H,{t))Gs{t,t') = 5{t~t')-s^ 1,2,3- --oo. (10) 

Eq. pl}|) can be solved in terms of the hnearly indepen- 
dent complete set of fundamental solutions, \(j)^j\t) >, 
for each s and all j, of the homogeneous Schrodinger 
equation, 



d 



{ih^^ - Hsit)Mf it) 0; s = 1, 2, 3, • • • , ^; allj. 



Explicitly, we have: 



(11) 



(12) 



allj 



We may now expand the total Green's function, 
GMcit,t')i using any basis Green' function, Gs(t,t'), in 
an infinite series, 

GMc{t,t') = Gs{t,t') + [ dhGs{t,h)DMc{ti,t') 



X Gs{ti,t') + / / dtidhGsitM) 



X DmC ih,t2)Gsih, t2)DMC {t2 ,t')Gs {t2 , t') 

+ •■•;s = 1,2,3,- ■ -,00. (13) 

where we have introduced the (inhomogeneous 
Schrodinger) operator: 

DMc{t, t') = [{HMc{t) - ih—) + iU{t - t')] (14) 

which, it is worth noting, depends on the total Hamil- 
tonian HMcit)i Sind not on its various partitions. The 
symbol J above stands for the integration over the entire 
time axis: / = J^^- We may note, parenthetically, that 
the presence of the theta-function in the explicit repre- 
sentation of a Green's function, Eq. p^ . automatically 
accounts for the appropriate domains of the intermediate 
time integrations in Eq. p^ . Eq. is a solution of 
the Green's equation ([8|). This can be easily verified as 
follows: First, we apply (fHl) on and use (O to get, 

DMcit, t')Gsit, t') 

- [{HMc{t) - ifi^) + iri5{t - t')]Gs{t, t') 
ot 

= [{H, (t) ~ iri^)Gs [t, t') + Vs {t)Gs {t, t') 
+ ih5{t-t')Gsit,t')] 

= [-d{t - t') + Vs{t)Gs{t, t') + 5{t - t')e{t ~ t') 

X Y.\c^^;\t)><c^^;\t')\] 

allj 

= {-bit - t!) + Vs (t)G(,) [t, t') + S{t - t')0{t - t')] 
= Vsit)G,it,t');s^l,2,3---,oo. (15) 



where, Vs{t) = HMc{t) — Hs{t), and in the last step 
we have used the fact that an integration variable, say 
tn+i, in any of the integrations in the series p3p . always 
approaches its upper limit, say as: lim(t„ — in+i) = 
0+, 9{tn — tn+i) — 1. Second, substituting Eq. p3|) in 
Eq. (I13|) we may sum the series on the right hand side 
as, 

GAic{t, t') 

= Gs{t,t') + J dtiGs{tM)Vs{ti) X [Gs{ti,t') 

+ / dt2Gs{tut2)Vs{t2)Gs{tlX) + ■■■] 



G,{t,t')+ / dtiGs{tMVs{ti) X [GMc{ti,t')] 



(16) 

where, we have used the fact that the quantity in the 
square brackets above is equal to the series itself. Next, 
operating on the last equation from the left with [iTi^i — 
Hs{t)), noting Eq. PH)) . and carrying out the resulting 
delta-function integration at once, we get, 

iih^^ - H,{t))GMc{t,t') - d{t - t') + Vs{t)GMc{t,t') 

(17) 

Finally, on transposing the last term to the left hand side 
and using Eq. ([9]) we arrive at, 

[in^^ - HMc{t))GMc{t,U) - 5{t - t') (18) 

which agrees with Eq. ([5]); q.e.d. Thus, Eq. p^ pro- 
vides a general series solution for the total Green's func- 
tion GMc{tj t') for an arbitrary choice of the basis Hamil- 
tonian Hs{t), and the associated basis Green's function 
Gs{t,t'). We may choose. 



Hs{t) = Hxit) 



(19) 

where A is a real number. The associated Green's func- 
tion G\(t,t') is defined by 

{ih^-Hxit))Gxit,t')^Sit-t') (20) 
ot 

We find its exact solution to be, 

Gx{t,t') = -j0(t-i')e'^'^^^*^" 



> e 



< 



allj 

= g^A^A(t).r^J^^^/)g-^A^A(^').r' (2I) 

The validity of the above solution can be established 
without difficulty by substituting Eq. ([2T|) in Eq. (|20|) . 
and simplifying by noting the definition ([7]) and using. 



(Pop - A-A(t))2e»^'fe^W- = e^^'fe^W"-p2 (22) 
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We now substitute Gs{t,t') = Gx{t,t') in Eq. ^ for 
the total Green's function, to get: 

GMc{t, t') 

= iGx{t,t')+ f dhG)^{t,t,)DMc{tut')Gx{h,t') 



dt2dhGx{t, h)DMc{ti,t2)Gx{ti,h) 
X DMc{t2,t')Gx{t2,t') + ---) 

X Gaitl,t')+ [ [ dhdt2Gait,h)D'^l^^{h,t2) 



X Ga{h,t2)D'^^^cit2,t')Ga{t2, t') + ■ ■ .Jg-A- A(t')-' 

(23) 



where, 



{[ 



(Pop + (A-l)fA(0) ^ 
2to 



■ + Va{r) + X-A{t)-r 



in—] + in6{t-t')} 



{[(^ + K(r)) - ^h-] + Vx{t) + zh5{t - t')) 

(24) 



and, 



Vxit) = [(A-1)— A(i)-Pop + (A-l)^-^A2(<) 



\-A{t)-r] 



(25) 



Operating with D\,jfj{t,t') from the left on to Ga{t,t'), 
noting Eq. ([5]), and calculating similarly as in the case 
of Eq. (fT5)) above, we get: 

= {{i^ + Va{r))-^h^^]Ga{t,t') 

+ Vx{t)Ga {t, t') + iM{t - t')Ga {t, t') } 

= Vx{t)Ga{t,t'). (26) 
Finally, using Eq. (^5)) in Eq. (|^ . we obtain 

GA/c(i,t') = e'^jfe^W-'-[Ga(i,<') + J dtiGa{t,h)Vxiti) 

X Ga{tut')+ [ I dhdt2Ga{t,h)Vx{tl) 

(27) 

This is a one-parameter family of infinite series represen- 
tations of the total Green's function GMci'tit')^ for a-ny 
value of the real parameter A : 



X Ga{tut2)Vxit2)Ga{t2,t') 



The transition amplitude, S'[^{tf,ti), between the 
field- free reference states 61: 



(28) 



and 

<^'fitf)\ = <</.^.(i/)|e-^^^^(*/)-'- (29) 

that are prepared initially at t' — ti and detected finally 
at a later time t = tj, where A{tij) are arbitrary con- 
stant vector potentials, is: 



= ih<(t>)(tf)\GMc{tf,U)\4>^{U)> 

= % dh < rfihWxihMiti) > 

+ J Jdhdt2 < 4>}{tl)\Vx{tl)Ga{tut2)\4>1{t2) > 

+ •••]; (A, a real number) (30) 

where, we have used the total Green's function G{tf,ti) 
given by Eq. ((27)) . cancelled the phase factors depending 
on the arbitrary (constant) vector potentials, and simpli- 
fied by putting, 

<rfmGaitf,h) = -l<c/,^(ti)| (31) 

and 

Ga{t^,UMiU)> - -^\^ntn)>. (32) 

Note that the above expression ([50)1 holds for any value 
of the real parameter A, and thus constitutes a "mas- 
ter" expansion of the transition amplitude derived in the 
minimal-coupling (MC) "velocity" gauge. It provides a 
one-parameter family of infinitely many series representa- 
tions or branches of the amplitude function Sl^(tf, ti). 
We may choose, for instance, the parameter A = 0, 1, or ^ 
on the right hand side of Eq. and get, respectively, 

the three expansions of the transition amplitude: 



= 5 



7i 



dii <(/.^(tl)|Kei.(tl)|C(^l) > 



dt^ < rf{tl)\V,eL{h)Gaiti,t2Mtit2) > 

;A = 

i 

" n 



(33) 



Sf^ " ^[ / dti < rfihWlenXtlMtitl) > 



+ J J dtidt2 < rf{tl)\Vlen.{h)Gaitl,t2) 

X Vlen.{t2Mf{t2) >+■■■]: X=l (34) 



dti < rf{t,)\Vhyb.{ti)\^'}{ti) > 
dtidt2 < rfiti)\Vhyb.{tl)Ga{tl,t2) 

X Vi,ybXt2Mtit2) > + ---];A = i 



(35) 
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where, from Eq. ([25)) . 

2 

Kei.W = (-— A(t)-pop + -^A2(t));A-0 

TOC Zmc 

VienXt) = (-A(i)-r) = (-e£;(i)-r);A = l 

c 

(36) 

The first two series (|33|) and (|34|) are readily recognized, 
on comparing the expressions of the respective interac- 
tions in Eqs. ((36)) . as the two well-known perturbation 
series, traditionally obtained in the "velocity" and the 
"length" gauge, respectively. Thus, they are seen to be 
nothing but two equivalent infinite series representations 
(branches) of the same amplitude function Sl^{tf,ti) 
and, hence, are but equal in their common domain of 
convergence (assuming, it exists) 0. We may add that 
at present, little, if any thing, is known regarding the con- 
vergence properties of these infinite series representations 
of the same amplitude function. The third expansion 
psp provides another equivalent series representation (in 



terms of a "hybrid interaction" ) of the same amplitude by 
choosing A = i, and so on for any other equivalent series. 
Before concluding, we point out that there is no diffi- 
culty in deriving a "master" expansion, analogous to Eq. 
([30)) . from the total Hamiltonian given in the "length" 
gauge (or for that matter in any other gauge) by pro- 
ceeding exactly analogously as shown above - it leads to 
the same conclusion of the mathematical equivalence of 
the two traditional perturbation series, in the velocity 
and length gauges, in a common domain of convergence, 
as demonstrated above. 

To conclude, starting explicitly in the minimal- 
coupling (MC) "velocity" gauge, we have derived a "mas- 
ter" perturbation expansion (PD|) for generating a one- 
parameter family of equivalent infinite series represen- 
tations of a given transition amplitude. The result is 
used to demonstrate (Eqs. p3|l and (p4| ) that the two 
well-known perturbation series, traditionally obtained in 
the "velocity" and "length" gauges, are only two equiva- 
lent infinite series representations or branches of the same 
amplitude function Sl^{tf,ti), and hence, are equal in 
their common domain of convergence (provided that the 
latter exists). 
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